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ABSTRACT 

We study the emission of X-rays from lobes of FR-II radio galaxies by inverse Comp- 
ton scattering of microwave background photons. Using a simple model that takes 
^^ , into account injection of relativistic electrons, their energy losses through adiabatic 

expansion, synchrotron and inverse Compton emission, and also the stopping of the 
jet after a certain time, we study the evolution of the total X-ray power, the surface 
,—, . brightness, angular size of the X-ray bright region and the X-ray photon index, as 

O^' functions of time and cocoon size, and compare the predictions with observations. We 

find that the radio power drops rapidly after the stopping of the jet, with a shorter 
£h ■ time-scale than the X-ray power. The X-ray spectrum initially hardens until the jet 

^ ' stops because the steepening of electron spectrum is mitigated by the injection of 

d . fresh particles, for electrons with 7 ^ 10 3 . This happens because of the concurrence 

of two times scales, that of the typical jet lifetimes and cooling due to inverse Comp- 
ton scattering (~ 10 7-8 yr), of electrons responsible for scattering CMB photons into 
£> ■ keV range photons (with 7 ~ -y/l keV/kTcMB ) • Another finding is that the ratio of 

Q\ ' the X-ray to radio power is a robust parameter that varies mostly with redshift and 

ambient density, but is weakly dependent on other parameters. We also determine 
the time-averaged ratio of X-ray to radio luminosities (at 1 keV and 151 MHz) and 
find that it scales with redshift as oc (1 + z) 38 , for typical values of parameters. We 
iy-\ . then estimate the X-ray luminosity function of FR-II radio galaxies and estimate the 

f^ ' number of these diffuse X-ray bright objects above a flux limit of ~ 3 x 1CP 16 erg 

O ! cm~ 2 s _1 to be ~ 25 deg~ 2 . 
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1 INTRODUCTION of the relativistic plasma in the lobe. Combined with ra- 

dio observations, this emission provides constraints on the 
A number of radio galaxies have been observed to emit dif- ,. n ,, , ,, ,. , ., ,. ,, , 

magnetic held and on the energy distribution among the rel- 

luse X-rays in the region between the nucleus and radio hot , . . ,. ,. , , ,, ,„., ., , . jn _ ,„4 . 

ativistic particles, or both. While electrons with 7 > 10 to 
spots in recent years. This emission has been interpreted as .. „ TT , ,. .. ,,. . , ,. „ , , 

J emit GHz synchrotron radiation (tor typical magnetic helds 

inverse-Compton scattering of the cosmic microwave back- c c „■. ., . lr .s ' ., „,,„ 

^ ,,..., or a few pG), it requires 7 ~ 10 to upscatter the CMB 



ground (CMB) photons by relativistic electrons in the radio 
lobes that have lost most of their energy and their Lorentz 
factor have come down to 7 ~ 10 3 (Fabian et al. 2003; Blun- 
dell et al. 2006; Erlund et al. 2008; Johnson et al. 2007; Cros- 
ton et al. 2005; Fabian et al. 2009; Isobe et al. 2009). That 
relativistic electrons in radio galaxy cocoons could upscatter 
CMB photons to X-rays has long been anticipated since the 
discovery of CMB (e.g., Felten & Rees 1967). 



photons to the observed 1-10 keV range in X-rays. 

The simultaneous radiation in two different wavelengths 
caused by the same parent electron population provides an 
excellent opportunity to probe the physical nature of the 
radiating object. Comparing the X-ray and radio emission 
from lobes, Blundell et al. (2006) showed that the injected 
electron energy distribution have a low-energy turn-over at 
7min ^ 10 3 , since the X-ray bright regions did not coincide 



Since the increase in the CMB energy density with red- ... ,. , .... . , 

4 . 1 with the radio hot spots. I he comparison between X-ray and 



shift (oc (1 + z) ) compensates for the cosmological surface 
brightness dimming, this extended emission has been seen 
at both low and high redshift, and has been used as a probe 



radio studies also shed light on the magnetic field. Croston et 
brightness dimming, this extended emission has been seen , ,„„„_% , ,, . ,, . r , . . „ , , , ,, . 

, , , ,,.,,,.. , , , , , al. (2005) argued that the interred magnetic held strength is 



close to the equipartition value (lying between 0.3-1.3 times 
the equipartition value). 

In this paper, we use a variant of a model of emission 
biman@rri.res.in from FRII radio galaxies advocated by Kaiser & Alexan- 
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der (1997) (hereafter KA97) and Kaiser, Dennett-Thorpe 
& Alexander (1997) (hereafter KDA97) to study the evolu- 
tion of X-ray power and surface brightness with time, and 
as functions of the radio lobe size, ambient density, redshift 
and the jet power. We do not assume self-similarity of the 
evolution of the radio lobe assumed by KA97 and KDA97, 
in order to study the effect of the cessation of jet activity on 
the X-ray surface brightness. 

The evolution of radio lobes after the stopping of jets 
has been discussed by Komissarov & Gubanov (1994), al- 
though they did not address the soft X-ray emission by old 
electrons. Kaiser & Cotter (2002) have also discussed the 
evolution of radio galaxies in connection with the observa- 
tions of radio relics, and Reynolds et al. ( 2002) have nu- 
merically studied the expansion of 'dead' radio galaxies. 

We first discuss the KA97 and KDA97 model and the 
changes introduced in the model in the present paper, and 
then discuss the results of our calculations. 



2 EVOLUTION OF RADIO GALAXY LOBES 
FOR FRII SOURCES 

2.1 The KDA97 model 

The KDA97 model of the dynamical expansion of an FR 
II- type radio source assumes a self-similar expansion, driven 
by twin jets emerging from the nucleus in opposite direc- 
tions, pushing the surrounding environment. The jets pro- 
duce strong shocks where the jet particles are accelerated 
and these particles made the cocoon expand. The density 
distribution in the ambient medium is assumed to be a 
power-law, p(r) — po{r/ao)~ , where po is the density at 
a core radius ao- If the jet power is denoted by Qj, then 
the combination (Qj/po%) 5 -P t 5 -f has the dimension of 
length, and, in a self-similar model, it is proportional to the 
length of the jet, Lj(t) (Falle 1991). Half of the cocoon vol- 
ume is approximated by a cylinder of length Lj and a base 
radius. The ratio R between the length and half-width of 
this cylinder is referred to as its axial ratio. 

The dynamics of the cocoon is determined by the co- 
coon pressure, p c — (T c — l)(u e + ub), where the contribu- 
tions from different components have been added: (a) rela- 
tivistic electrons with energy density u e and adiabatic index 
r c and (b) a tangled magnetic field with energy density ub- 
(This implicitly assumes that T c — 4/3, the adiabatic index 
appropriate for magnetic field pressure.) Suppose the cocoon 
volume scales as t ai . Then the cocoon pressure evolves with 
time as p c oc t~ r " ai (see KDA97). The cocoon is assumed 
to be divided into small volume segments filled with magne- 
tized plasma and particles that are injected into the cocoon 
at the jet terminal shock. 

The initial electron energy distribution is assumed to 
be a power law in energy, 



n('ji,ti)d'y = n j i p dji 



(1) 



where the spectrum ranges between 'y m in and ^ ma x- The 
electrons lose energy through adiabatic expansion, syn- 
chrotron and inverse Compton losses. KDA97 found a 
closed-form solution for the energy distribution at a later 
time n(7, t) in the case of a self-similar expansion of the 
cocoon. 



The ratio of the energy density in particles to that in 
magnetic field is assumed to be r, so that the magnetic en- 
ergy density ub at time t is given by, 

rpc(t) 



u B {t) 



(2) 



(T c -l)(r + l) 

Assuming that for synchrotron radiation, electrons emit only 
at their critical frequency v = 7 jl, where vl is the Lar- 
mor frequency, the total radio power emitted by electrons 
in a volume element SV of the cocoon (radio power per unit 
frequency and solid angle) is given by, 

3 

(3) 



1 
5P V — —(jtcub — nh)5V . 

07T V 



Here the volume element at time t is related to the initial 
time interval over which the particles residing in it were 
injected, as, 



5V(t) = 



(rc ' 1)Qj (4fi 2 ) (1 - r - )/r - 



8U. 



(4) 



The total radio power is then calculated by integrating over 
all time, 



(5) 



P v = / 8P V , 

Jo 

using eqns © and ((4]). 



2.2 Evolution without self-similarity 

In this paper, we would like to study the evolution of radio 
and X-ray power on a time scale that is longer than the 
typical jet lifetime, and study the effect of the stopping of 
jets on the X-ray emission. As we will explain later in the 
paper, the stopping of the jet has important effects on X-ray 
power of the cocoon, and the results are markedly different 
than in the case where the jet continues to be active. To 
do this, we cannot use the self-similar evolution assumed in 
KDA97. Instead, we calculate the size evolution the radio 
lobe in the following way that captures the basic physical 
processes in a simple manner. 

Following the standard evolutionary picture (Scheuer 
1974; Begelman & Cioffi 1989, Reynolds & Begelman 1997), 
that radio jets inject relativistic particles into a cocoon, 
which is overpressured (Nath 1995) and drive a strong 
shock into the ambient medium. The speed with which this 
shocked shell expands is determined by the ram pressure of 
the ambient gas that enters the shock. Assuming that the 
pressure inside the cocoon at any given time is uniform (see 
KA97), and a value p c (t), and assuming that radiation loss 
is small, we have for the expansion of the cocoon, 

Qj(t) = -f T{V c p c + Y c p c V c ) 



dt 



T c - I 

PO 



(6) 



where V c is the volume of the cocoon, and po is the ambient 
density. For simplicity, we assume a constant axial ratio R, 
and the volume is assumed to be, V c = t§t£j (KA97). 

For the jet luminosity we specify a jet lifetime tj, so 
that Qj is a constant for t ^ tj, and Qj = for t > tj, 
where tj ~ 10 7 8 yr. Numerical simulations of 'dead' radio 
galaxies after the cessation of jet activity have shown that 
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the overpressured cocoon continues to expand until a pres- 
sure equilibrium is established with the ambient medium 
(Reynolds, Heinz, Begelman 2002). 

We assume a constant axial ratio R ~ 2, which is an av- 
erage value (Leahy and Williams 1984), and, following Wang 
& Kaiser (2008) , we use a profile of the ambient density as 
given by, 



Pa = po(r/ao)~ 



Ar 



(7) 



where A = poa^. For ft = 2, A has the units of [g cm -1 ], and 
one can infer its value from the observations of environments 
of radio galaxies. Typically, for ft = 2 and ao ~ 200 kpc, an 
electron density n e ,o(= Po/(pm p ) ~ 10" 4 cm -3 would imply 
A ~ 10 19 g cm" 1 , a value we assume for most of our models 
in this paper. Observations of individual galaxies show A to 
be in the range of 10 19_2 ° g cm -1 (Fukazawa et al. 2004), 
whereas Jetha et al. (2007) inferred A ~ 10 19 -10 20 g cm" 1 in 
group environments, although these determinations assume 
ft ~ 1.5. For simplicity, because of the ability of combining 
two free parameters (po and ao) into one (A), we assume a 
value of ft — 2 for our calculations. 

Solving equations [6l we can therefore determine the co- 
coon pressure p c (t) at any instant, from which we determine 
an equipartition value of magnetic field, using eqn ^. as- 
suming (as in KDA97) a ratio r between particle and mag- 
netic field energy density. The corresponding electron en- 
ergy density is then given by u e (t) = u B (t)/r, assuming 
that there is no thermal particles in the cocoon (i.e., in the 
language of KDA97, ut = 0). 

In the KDA97 model, particles are assumed to be in- 
jected at the termination point of the jet. In our case, after 
the jet stops, the injection of particles also stops, and so the 
integration in the equation (0, for the determination of ra- 
dio power at a given time t, has an upper limit of min(i, tj). 
To perform this integration, we rewrite the equation for vol- 
ume segments (eqn @) which in the KDA97 model used the 
exponent oi for self-similar evolution of cocoons. For t ^ tj, 
any small volume segment of the cocoon at a a given instant 
t can be related to the pressure at the time (ti) when the 
electrons in this segment were injected into the cocoon. One 
can rewrite equation Q in the KDA97 model as (for t < tj) 



8V(t) = 



(r e -l)Qj ( 2 )( i-r e )/r e 



( Pc(t) \ 
\Pc(ti)J 



i/r c 



5U. 



(8) 



The evolution of Lorentz factor is explicitly solved using 
the loss equation, 



d7 11 dV c 4 a T i, 

17 = ~o77~rr ~ q — T ( Us +Uc 
at 6 Vc at A m e c 



0) 



where the first term denotes adiabatic energy loss and is 
computed using the results of equations ©). The second 
term combines radiation loss in synchrotron and inverse- 
Compton scattering. Here u c — oTq MB is the CMB photon 
energy density, or is the Thomson cross-section, m e is the 
electron mass and c is the speed of light. 

Thus, beginning with an initial energy distribution law 
with a power-law index p, one can solve for the energy dis- 
tribution at any given time, 71(7, t), given the initial dis- 
tribution, 71(74, ti). Note that this was analytically done by 
KDA97 for a self-similar evolution of the cocoon, and we 
explicitly solve it in order to go beyond self-similarity. 

Using this knowledge of 71(7, t), we can then use equa- 



tion ((3]) in conjunction with eqn (J5J) , and integrate over time 
to calculate the radio power at a given frequency and at a 
given instant. 



2.3 Inverse-Compton radiation 

We extend this formalism further to determine the inverse- 
Compton (IC) emission as a function of time. Firstly, we 
note that electrons with Lorentz factor 7 boosts a CMB 
photon of frequency vcmb into an energy h^vcMB- The 
precise calculation of the IC power for this photon would 
require one to consider the total spectrum of CMB pho- 
tons. But we can simplify it for our purpose here by as- 
suming all CMB photons to have a single frequency. We 
assume that the CMB photon distribution function is given 
by v'(e) = v' 5(e — knTcMB), where the normalizing factor 
v'o can be calculated by requiring that total energy density 
f ev'{e)de = aT^ MB , and one has v' — aTQ MB /kB- The 
total scattered power depends on the integral (see equation 
(7.29a) in Rybicki & Lightman 1979), 



I' 



dev (e) e" 



an. 



kn 



'-{kaToi 



\(p-i)/2 



(10) 



If one had used the blackbody distribution function v(e) = 
(87re 2 //i 3 c 3 )(exp(e/fcsToAfs) — l)" 1 , then the corresponding 
integral would have yielded (see equation 7.31 in Rybicki and 
Lightman 1979), 



dev{e) — = JL_ {kBTcMB) — T{ P±l )a 



P + 5, 



,(11) 



where the symbols have the standard meanings. Then the 
ratio, 



1 



15 r(*±5) C(2±5) ' 



(12) 



shows the error one incurs in assuming all CMB photons to 
have the peak frequency in estimating the inverse Compton 
power. For p — 2.2, one finds I 1 / 1 ~ 1.6. Therefore the total 
scattered power calculated using this assumption is correct 
within an accuracy of 60%, and we adopt it for simplicity in 
our calculation. 

We therefore calculate the IC power at a given fre- 
quency, in the manner of previous eqn ([3]), as, 



Pic, v 



i[t,*il 



1 y 

arcuc — n{^)8V , 



6tt 



(13) 



using eqn ((8| for SV , and the value of 71(7,4)) found from 
solving eqn @. 

Not all volume elements in the cocoon would contribute 
to the IC radiated power though. Electrons in some volume 
element that was injected at an earlier time ti may lose 
energy to the extent that they fail to boost CMB photons 
to the observed X-ray band. These volume elements would 
not contribute to the X-ray power. 

We calculate the total projected area of the cocoons 
that contribute to X-ray power, assuming the cocoon to be 
in the plane of the sky, and an axial ratio R. This allows us 
to calculate the X-ray surface brightness. We also calculate 
the fraction f x of the total surface area of the cocoon that 
contributes to X-ray emission. 
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Figure 1. Time evolution of radio and X-ray powers are shown 



in this plot. Thick lines show the result for Case I (Qj 



10 45 



thin lines show the results for Case II ( 



10 erg s 



z = 1, with the same ambient density as in Case I). Dotted, solid 
and dashed lines show the cases for t 3 ■ = 5 X 10 7 , 10 8 and 10 9 yr, 
respectively. Time evolution of the cocoon size (Lj ) is shown in 
the top-left panel, radio power at 178 MHz (in W/Hz sr) in the 
top-right, and the inverse Compton power at 1 keV in the same 
units is shown in the bottom -left panel. The bottom right panel 
shows the corresponding ratio of inverse Compton to radio power 
as a function of cocoon size for the three cases. Two horizontal 
dot-dashed lines indicate the range of observed ratios (see text). 



3 RESULTS 



To begin with, as a fiducial case, we consider Qj = 10 
erg s _1 , A = 10 19 g cm -1 at z = 0.1, and a constant axial 
ratio R — 2. Following KDA97, we also adopt r — 0.785 and 
p — 2.14, "fmax = 10 6 and 7 m i„ = 1. We refer to this set of 
parameters as Case I. 

The set of thick lines in Figure 1 show the results for 
this fiducial case. We plot the time evolution of cocoon size 
in the top-left panel, radio power at 178 MHz in the top- 
right panel, IC power at 1 keV in the bottom-left panel. The 
corresponding evolution of the ratio of synchrotron to IC 
power with cocoon size is shown in the bottom-right panel. 
The solid lines show the case for tj — 10 8 yr, and dotted 
and dashed lines show the cases tj = 5 x 10 7 and 10 9 yr, 
respectively. 

In the same figure, we also show the results of another 
case, with a higher jet power and one located at a higher 
redshift than in the fiducial case. The thin set of lines in 



Figure 1 show the case with Qj 



10 



16 erg s 1 at z 



1 



(Case II), keeping the values of other parameters the same. 
The evolution of the size of the cocoon is self-similar, 



with sources at higher redshift dimming faster than their 
low redshift counterparts because of inverse Compton loss. 
For shorter jet lifetime, the radio luminosity drops precipi- 
tously after t ~ tj. This is mainly driven by rapidly declining 
pressure inside the cocoon owing to adiabatic expansion. 

The X-ray luminosity of the cocoon at 1 keV (bottom- 
left panel of Figure 1) increases till t ~ tj, as high energy 
electrons lose energy (dE/dt oc E 2 ) and enter the energy 
bins that are conducive for inverse Compton scattering CMB 
photons to the soft X-ray band. The number of electrons 
with Lorentz factors ~ 10 3 increase because the draining of 
these electrons to lower energy bins is more than compen- 
sated by the supply of energetic electrons from the jet. But 
after the jet stops, no new electrons are supplied to the co- 
coon, and the X-ray power drops rapidly (more so at high 
redshifts). 

The time scale for inverse Compton loss is tic ~ 3 x 10 7 
yr for 7 ~ 10 3 at present epoch. The time scale for the drop 
in X-ray power is made shorter than this by the additional 
loss of energy from adiabatic expansion of the cocoon. It is 
clear from curves in Figure 1 that the stopping of the jet ac- 
tivity has a profound effect on the X-ray power that would 
have been missed in a calculation with self-similar evolution 
of cocoons. Self-similar models would have predicted a ris- 
ing X-ray power, although the rate of increase would have 
tapered beyond a time-scale ~ tie, because the loss of elec- 
tron energy would have been continually compensated by 
injection of new and energetic electrons. In the case of jet 
stopping after tj, this privilege is withdrawn from the co- 
coons and their X-ray power rapidly diminishes with time. 

The ratio of X-ray to radio power (bottom-right panel), 
therefore, shows a mild increase till tj, after which it in- 
creases rapidly, because of the drop in radio luminosities. 
We note that, typically the value of the ratio at t ~ tj is 
of order Pi C /P s ~ 10~ 10 -10 -9 . For a single power law elec- 
tron energy distribution with index p — 2a + 1, where a is 
the synchrotron (and inverse Compton) radiation spectral 
index, the ratio of relative luminosities can be written in 
terms of the strength of the magnetic field and the CMB 
temperature. For p = 2.5, we have (e.g. Tucker 1977), 



p. 

1 IC 



4.7 x 10" 



11 ( TcMB 

2.73 K 



3+ct 



B 



10 nG 



-(<*+!) 



(14) 



where the coefficient is ~ 2.7 x 10 -11 for p — 2, and 
8.3 x 10 _1 for p — 3. These values are somewhat lower than 
predicted from our detailed calculations, which take into 
account the electron energy distribution in different parts 
of the cocoon as they evolve. Although the magnetic field 
strength throughout the cocoon is assumed to the same at 
any given time, the difference in the evolution of electron 
energy spectrum in different parts of the cocoon can give 
rise to a value of the ratio P x / ' P r that is different from the 
above simple estimate. Admittedly, the flux ratio is mainly 
driven by magnetic field, which is what the above simple 
estimate indicates, but our calculations show the important 
role played by the evolution of electron spectrum prior to 
the observed epoch for an object. 

The observed ratios of fluxes at 1 kev and 178 MHz fall 



as expected, till t ~ tj (top left panel of Figure 1), after in the range of P x /P r 



10- 



-10" 



5 , for objects with 



which the cocoon size expands slower than before. The radio 
luminosity as a function of time (top-right panel) is compa- 



rable to the results of KDA97 for the case of tj 



10 a 



yr, 



a variety of linear sizes and at different redshifts (e.g., 3C 
215 (2 = 0.41), 3C 334 (z = 0.56), 3C 9 (z = 2.0) (Bri- 
dle et al. 1994), 3C 219 (z = 0.17) (Perley et al. 1980), 
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3C 265 (z = 0.81) (Bondi et al. 2004), 3C 109 (z = 0.31) 
(Giovannini et al. 1994), 3C 179 (z = 0.85) (Sambruna ct 
al. 2002)). We indicate this range by two dot-dashed lines in 
the figure. Since the evolution of this ratio depends on many 
parameters (as listed in Table 1, for example), our goal in 
this paper is not to reproduce all the observed parameters, 
but to probe whether or not the model yields representative 
values consistent with observations. 

It is apparent from the figure that the observed range 
of ratios are mostly consistent with model predictions for 
different parameters. It also appears that the observed range 
is not related to the rapidly rising part of the evolutionary 
curves. In other words, the systems observed in X-ray (near 
1 keV) should also be radio bright at low frequencies, of 
order ~ 150 MHz, that we have used for our calculations 
(~ 178 MHz). Low frequency radio observations therefore 
should be an important probe of X-ray bright cocoons, as 
has been pointed out by previous authors. 

Another aspect is that cocoons spend a relatively short 
time in the rapidly rising parts of the curves for ratio be- 
tween Px/Pr, because of the precipitous drop in radio lu- 
minosity after tj. In other words, radio galaxies that have 
faded from the window of frequencies of order a few hundred 
MHz, are unlikely to be X-ray bright, since the X-ray power 
would also drop (see bottom- left panel of Figure 1). 

Next, we plot in Figure 2 the corresponding flux in the 
1-5 keV band (top left panel) for the cases shown in Figure 
1. Our results show that for the parameters used here the 
cocoons typically emit at a flux of ~ 10~ 14 erg s— 1 cm~ 2 
per keV (or a few micro-Crab) before dropping after tj , con- 
sistent with observed fluxes (see, e.g., Laskar et al. 2010). 

We also need to consider the X-ray surface brightness 
apart from the total flux. We calculate the total area of 
the X-ray bright region of the cocoon, by summing over 
the volume elements in which electrons contribute to the 
X-ray emission in the 1-5 keV band, and projecting in the 
plane of the sky. This is shown in the bottom-left panel of 
Figure 2, which suggests that for cocoons larger than 100 
kpc, typically a patch as large as tens of arcseconds would 
be X-ray bright. If one estimates the fraction of the cocoon 
area that is X-ray bright, this fraction initially rises (till tj) 
to ~ 0.3 — 0.5 of the total area (assuming the cocoon the 
projected in the plane of the sky), and then drops. 

This calculation of the angular size of the X-ray bright 
region of the cocoon allows us to estimate an average X-ray 
surface brightness of the cocoon, and we show the results in 
the top-right panel of Figure 2, for the same cases as men- 
tioned earlier. The surface brightness (in the 1-5 keV band) 
initially drops in a gradual manner, owing to two competing 
effects: increasing X-ray luminosity and an increasing por- 
tion of the cocoons which are illuminated by X-ray. But the 
it drops rapidly after tj, especially at high redshift. 

Another important probe of X-ray bright cocoons is 
the X-ray spectral index. We calculate the photon index 
r(= a + 1) in the 1-5 keV band for the same cases, and the 
results are shown in the bottom-right panel of Figure 2. We 
also include data points for a number of cases which have 
been observed long enough for spectrum determination. We 
have estimated the cocoon angular size from visual inspec- 
tion from the observations cited in the caption, assuming 
the cocoons to be projected in the plane of the sky, and 
estimated the physical size assuming a A-CDM cosmology 
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Figure 2. Evolutions of X-ray flux in the 1-5 keV band (top-left 
panel), X-ray surface brightness (top-right), angular size of the 
X-ray bright region (bottom-left) are shown as functions of the 
cocoon length, for the same cases as in Figure 1. Thick lines refer 



to the fiducial case of Qj = 10 45 erg s 



A 



z = 0.1, with 7„ 



1 (Case I). Thin lines consider 



-i 
10 46 

erg s — 1 at 2 = 1 (Case II). Solid, dotted and dashed lines show 
the cases for t 3 • = 5 X 10 7 , xlO 8 and 10 9 yr, respectively. The 
bottom right panel shows the corresponding X-ray photon index 
in the 1-5 keV band as a function of cocoon size for these cases. 
We also plot the data points for 3C 47N (z = 0.43, Lj ~ 333 
kpc), 3C 215 (z = 0.41, Lj ~ 246 kpc) (Bridle et al. 1994), 3C 
219 (z = 0.17, Lj ~ 270 kpc), 3C 452 (z = 0.81, Lj ~ 45 kpc) 
(Perley ct al. 1980), 3C 265E (z = 0.81, Lj ~ 55 kpc) (Bondi ct 
al. 2004), for comparison. 



with h = 0.7, Oo = 0.3, Qa — 0.7. The error bars on the 
data points are large, but there appears to be a trend of 
decreasing photon index with increasing cocoon size: larger 
cocoons show flatter X-ray spectrum. 

Our model predictions show that initially the X-ray 
spectrum gradually gets flatter (smaller value of F) till 
t ~ tj. This is because of the fact that for electrons 7 ^ 10 3 , 
the injection of fresh particles mitigates the radiative steep- 
ening of electron spectrum. The injection of electrons from 
the jet compensates for synchrotron energy losses for elec- 
trons with 7 ~ 10 3 , where two relevant time scales become 
comparable: (a) the jet lifetime (10 7 8 yr) and the (b) in- 
verse Compton cooling time scale of electrons with 7 ~ 10 3 
(~ 3 x 10 7 yr at present cosmological epoch). This energy 
scale (corresponding to 7 ~ VlkeV UTcmb ~2x 10 3 ) also 
happens to be the one required to upscatter CMB photons 
to keV range. The concurrence of these two timescales pro- 
duces a temporary hardening of the soft X-ray spectrum. 

Then, after the jet stops, the number of electrons at this 
energy scale rapidly decreases because of inverse Compton 
and adiabatic loss, and the spectrum becomes soft again. 
This softening occurs rapidly at high redshift, as expected 
from increasing inverse Compton loss. With regard to some 
of the data points with large cocoons and T ~ 1.5, we found 
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Model 


Qj (erg s" 


- 1 ) 


2 


tj (yr) 


A 


(g cm 


- 1 ) 


7mm 


Case I 


10 45 




0.1 


10 8 




10 19 




l 


Case II 


10 46 




0.1 


10 8 




10 19 




l 


Case III 


10 46 




0.2 


10 8 




10 19 




l 



Table 1. The values of different parameters used in the models 
cited in the text are tabulated here. Apart from these parameters, 
all models use /3 = 2, T c = 4/3, R = 2. 



that they can be explained with models using a jet power of 
order 10 46 erg s _1 at low redshifts. 



3.1 Variations with parameters 

Although we have shown the results for cocoons for a few 
cases, varying parameters such as the jet power, redshift 
and jet life time, there are other free parameters in this 
model whose effects must be understood. The X-ray prop- 
erties of cocoon also depend on the ambient density (here 
parametrized as A) and the lower cutoff in the electron en- 
ergy distribution (j m in)- To study the effect of these pa- 
rameters, we first plot in Figure 3 (with solid lines) the 
results for a fiducial case (Case III): Qj = 10 46 erg s _1 , 
z = 0.2, tj = 10 8 yr, and A = 10 19 g cm" 1 , ^ mm = 1. In 
this section we will study the effect of changing the values of 
Lj, A, 7 m i„ separately, and compare with the results for our 
fiducial case. We will not vary the redshift and the jet life 
time, as the discussion in the previous section has already 
considered the effects of changing these parameters. Also, 
we keep ^ max fixed at 10 6 , and we do not change the values 
of p and r. 

In Figure 3, we have plotted the X-ray flux (in 1-5 
keV band; top-left panel), ratio between X-ray and radio 
power (top-right panel), X-ray surface brightness (in 1-5 
keV; bottom-left panel), and the photon index (r) in the 
bottom-right panel. All these parameters are plotted against 
the cocoon size. 

First, the results for the fiducial case mentioned above 
are shown with solid lines in all panels for easy comparison. 
Then, with dotted line, we show the results of changing jet 
power to 10 47 erg s _1 . As expected this cocoon expands to 
occupy a large volume, and is more X-ray bright than the 
fiducial case. The ratio of X-ray to radio power does not 
change, however, but the surface brightness increases. The 
photon index continues to decrease and the X-ray spectrum 
becomes harder with time (and size) for t ~ tj after which 
it becomes softer. 



Next, we change the value of r y m i n to 10 3 (keeping all 
other values same as in the fiducial case), and the results 
are shown with short-dashed lines. Increasing the lower cut- 
off results in increasing the number density of electrons for 
a given total energy in the cocoon, and produces a larger 
X-ray flux and surface brightness. But since there are more 
electrons with high energy in this case compared to the fidu- 
cial case, they also lose energy faster, and consequently, the 
drop in X-ray power after the stopping of the jet is faster 
than in the fiducial case. Interestingly, the increasing of 7 m i„ 
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Figure 3. Evolutions of the X-ray flux (in 1-5 keV band; top- 
left panel), ratio of X-ray to radio power (top-right), X-ray surface 
brightness (bottom-left) and X-ray photon index (bottom-right) 
are shown as functions of the cocoon size for a few cases. The solid 
line shows the results for Case III: Qj = 10 46 erg s _1 , A = 10 19 g 
cm -1 at z = 0.2, with 7 m i n = 1, and tj = 10 8 yr. The dotted line 
shows the case for Qj = 10 47 erg s _1 , keeping other parameters 
fixed. The short-dashed lines refer to the case with 7 m ; n = 10 3 , 
and the long-dashed lines show the effect of changing the ambient 
density to A = 5 X 10 19 g cm -1 . Data points from Figure 2 for 
photon index are again plotted here. 



does not effect the the X-ray to radio power ratio and the 
X-ray photon index. (It would, if the lower cutoff were to 
increase beyond 7 ~ 10 3 .) 

Finally, we change the ambient density to A = 5 x 10 19 g 
cm -1 that is indicative of a dense environment. The results 
are shown with long-dashed lines, and they show that the 
radio cocoon is constrained to remain small in size in this 
case. This results in not allowing the cocoon pressure (and 
magnetic field) to decrease fast enough, so that electrons 
are made to lose energy faster through synchrotron radia- 
tion. Both X-ray and radio power drop rapidly in this case 
and the ratio between X-ray to radio power also decreases 
(although, since the curves are drawn as functions of D, the 
rapid decrease is not apparent from the figures). The pho- 
ton index also hardens faster than in the fiducial case. This 
implies that it would be difficult to detect X-ray bright co- 
coons in dense environments like galaxy clusters, although 
it would be difficult to observe them anyway because of the 
bright X-ray emission from the intracluster gas. 

The X-ray fluxes shown in the top-left panel of Figure 
3 suggest that for typical radio galaxy parameters as used 
here, the cocoons shine in soft X-rays with a flux of 1-10 
micro-Crab, or, 10~ 14 -10~ 13 erg cm~ 2 s _1 per keV, before 
dropping to lower values after tj. 

Curves in the top-right panel of Figure 3 show that the 
ratio of X-ray to radio power is a robust parameter, in the 
sense that it does not vary with changes in jmin and jet 
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power at fixed times. The curves are plotted againstD, but 
the kink in the curves due to the stopping of jet occurs for 
all curves at the same time, namely, tj, and a comparison 
of the solid, short dashed (which happens to be superposed 
on the solid curve) and the dotted curves show that the 
ratio is a robust parameter. It does vary with redshift (as 
we have seen in Figure 1) and with ambient density (long 
dashed line in Figure 3). This robustness allows us to use 
this parameter to characterize the connection between radio 
and X-ray properties of cocoons, and we will discuss this 
issue further in the next section. 



3.2 Time-averaged ratio of X-ray to radio power 

Our results show that many of the X-ray properties of FR- 
II radio galaxy cocoons vary substantially with time, even 
when all other parameters such as jet power, ambient density 
and others are kept constant. It is therefore not easy to 
predict the X-ray properties of these sources that can be 
tested with observations, because it is difficult to determine 
the age of these sources from radio or other observations. 

One can however use the fact that the evolutionary 
time-scale of radio galaxies, of order a few hundred Myr, 
is much shorter than the Hubble time, even at redshifts 
when the radio galaxy population peaked in number density 
(z ~ 2). This implies that one can use time-averaged quan- 
tities related to X-ray emission, and speculate upon their 
average properties that could be observed and tested. 

In particular, we would like to determine the aver- 
age property of the ratio of luminosities in X-ray and ra- 
dio frequencies, v x P x /v T P r of the sources under considera- 
tion, averaged over time until their X-ray or radio emission 
drops rapidly: J ^j^ dtj J dt. Our results (Figure 3, top- 
right panel) shows that this ratio of luminosities is a robust 
parameter, and it varies negligibly with the variations in jet 
power, lower cut-off in 7, although it varies strongly with 
changes in ambient density. This ratio is however likely to 
increase with redshift, and we wish to determine the scaling 
with redshift. 

We therefore compute the X-ray luminosity v x P x in the 
1-5 keV band, and we choose a radio frequency of v r = 
151 MHz to compute v r P r . We compute a time-average of 
this ratio, summing over the duration in which the sources 
are both X-ray and radio bright in these frequencies (in the 
source frame). 

Celotti & Fabian (2004) have discussed the significance 
of this ratio of X-ray and radio luminosities in the case of 
radiation from electrons with a single power law energy dis- 
tribution. In this case, this ratio can be written as, 



Vx*x 



UCMB 

U B 



/ VxVB y 

\V r VCMB) 



(1 + *) 



3 + a 



(15) 



where a is the radio spectral index, Ucmb and Ub are the 
energy density in CMB photons (at z — 0) and magnetic 
field, respectively. Also, vb and ucmb are the gyrofrequency 
and the peak frequency of CMB photons at present epoch. 
For a field strength of B = 0.1-10 /xG, the ratio between the 
luminosities at 1 keV and 151 MHz (in the observer frame) 
can take values 0.08-300 at z = 0, for p = 2.6, and 0.08-50 
for p — 2. The ratio increases with redshift, although the 
increase can be mitigated by the possibility that magnetic 
field may change with redshift. Celotti & Fabian chose a 




Figure 4. Time-averaged ratio between X-ray and radio lumi- 
nosities is shown as a function of source redshift, for p = 2.5 
(solid line), and p = 2.14 (dotted line). The fit (~ 0.14(1 + z) 3 ' 8 ) 
is shown with a dashed line. 



conservative value of unity for this ratio, at all redshifts, 
and discussed the possible implications. 

In our model that includes evolving electron popula- 
tions, there is no single power law energy distribution of elec- 
trons, and the magnetic field is estimated through equiparti- 
tion arguments. We show the evolution of the time-averaged 
ratio of X-ray and radio luminosities at the above mentioned 
frequencies (in the source frame) with redshift in Figure 4, 
for two values of p — 2.5 (solid line), and p — 2.14 (dotted 
line), the last choice of p being motivated by KDA97. Inter- 
estingly, the difference between these two initial values of p 
is not large, and makes the time-averaged ratio of luminosi- 
ties a robust quantity. We find that the case for p — 2.5 can 
be fitted with a simple scaling with redshift, 



(i> x P x /v r P r ) t ~0.U(l + z ) 



(16) 



The results of the two cases of different values of p are differ- 
ent from the expectation from the simple formula mentioned 
above. The X-ray to radio luminosity increases with decreas- 
ing p (or increasing a), especially at high redshift. This is 
because the increased population of high energy electrons 
(for smaller values of p) rapidly lose energy, decreasing 7 
down to the level where they become X-ray bright, and this 
process become more efficient at high redshift. 



3.3 X-ray luminosity function of FR-II radio 
galaxies 

We are now in a position to estimate the number density of 
X-ray bright radio galaxy cocoons using the radio luminos- 
ity function of these sources, and using the above results for 
the relation between X-ray and radio emission. We use the 
radio luminosity function of FR-II galaxies, as determined 
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by Willott et al. (2001) from 3CRR, 6CE and 7CRS sam- 
ples. This luminosity function was determined assuming a 
cosmological model with fio = = $1a, fife = 1, h = 0.5. 
We have converted it for the A-CDM cosmology with £7a = 
0.7, fi m = 0.3, h = 0.7 using the relation (Peacock 1985), 



m dVi m -,dV 2 

Pl(Pl,z)— = P2(P2,Z)—- 

dz dz 



(17) 



where P is the luminosity derived in a specific cosmological 
model for a measured flux and at a given redshift z, and the 
indices refer to the two different cosmological models. The 
luminosities in two different models are related as, 



Pi£>r = Pi 

where 
*<*> = £ 



dz' 



y/n M {i + z')3 + n A 



(18) 



(19) 



is the comoving distance in a flat cosmological model. 
Also, the comoving volume in eqn [T7] is given by dV = 
4n D{z) 2 dD(z). 

Having done this, we have associated a diffuse X-ray lu- 
minosity to each FR-II sources with the scaling result from 
the last section, and determined the X-ray luminosity func- 
tion of FR-II galaxies. 

Figure 5 shows the computed X-ray luminosity function, 
in the units of number per unit (comoving) Mpc 3 , per 10 44 
erg s _1 , for three redshifts: z — (solid), z = 1 (dotted) and 
z — 2 (short-dashed line). The peak of the X-ray luminosity 
function shifts to increasing luminosity because of the strong 
redshift evolution of the X-ray to radio power ratio. We note 
that Celotti & Fabian (2004) assumed a constant ratio of 
unity and so their X-ray luminosity function peaked at the 
same luminosity at different redshifts (their Figure 3). 

We compare the computed luminosity function with 
that of clusters in X-ray in the 0.5-2 keV band (shown here 
with long-dashed line) as determined from the REFLEX 
sample by Bohringer et al. (2002) ( for a similar cosmol- 
ogy but with h = 0.5). Recently, Mullis et al. (2004) found 
that the cluster X-ray luminosity function decreases with 
increasing redshift at the high luminosity end, but does not 
evolve significantly below 10 44 erg s _1 . The comparison with 
the luminosity function expected from FR-II galaxies shows 
that the number density of these sources become comparable 
at the high luminosities at z ~ 2. 

We can also estimate the number of diffuse X-ray emit- 
ting cocoons in a given area of the sky above a given X- 
ray flux limit. Integrating the above luminosity function to 
z ~ 2, we find that the numer of sources above a flux limit 



gree. Recently Finoguenov et al. (2010) detected 
emitting radio lobes in 1.3 deg 2 above 2 x 10" 



X-ray 
15 erg cm -2 
and given the uncertainties, our estimate is consistent 



with it. 



4 DISCUSSION 

The expected number of X-ray bright radio galaxy cocoons 
as estimated above can be compared with the observed val- 
ues. Bauer et al. (2002) found six extended sources in the 
Chandra Deep Field North survey, within an area of ~ 130 
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Figure 5. Predicted X-ray luminosity function of FR-II radio 
galaxy cocoons, based on the relation between X-ray and radio 
power, is shown here, for z = (solid), z = 1 (dotted) and 2 = 2 
(short-dashed line). The X-ray luminosity refers to that in the 
1-5 keV band, and the luminosity function shows the comoving 
number density of objects, per 10 44 erg s — 1 , in the A-CDM cos- 
mology, with h = 0.7. The long-dashed line shows the present 
day X-ray luminosity function of clusters, based on the REFLEX 
(Rosat-ESO Flux Limited) survey by Bohringer et al. (2002). 



\ll deg 2 , at 
\ Our 



arcmin 2 , implying a surface density of ~ 167 

a limiting soft X-ray flux of ~ 3 x 10 -16 erg ci 

estimate shows that a small fraction of order ~ 0.1-0.3 of 

such extended soft-Xray sources in the sky could be due to 

FR-II radio galaxies. 

The X-ray luminosity function in the soft band can 
also be used to estimate mirror effect of X-ray emission, 
namely, the Sunyaev-Zeldovich (SZ) effect on the CMB. To 
some extent, it would underestimate the effect because of 
the fact the soft-band X-ray power is smaller than the to- 
tal X-ray power. Keeping this in mind, we can determine 
the ratio of the total radiation energy density in the soft- 
Xray band that is emitted by these sources and that is 
present in the CMB, by integrating the luminosity func- 
tion: ^ ~ l"*^ JL x dp ^' z) dL x . We find that 

Ae/e ~ 2 x 10~ 7 , integrating up to a redshift z ~ 2, al- 
though strictly speaking this is a lower limit. This implies 
a Compton y-parameter of order y ~ (1/4)— ^ 10 -7 . Ya- 
mada, Sugiyama & Silk (1999) considered the distortion of 
the CMB from the population of radio galaxy cocoons, us- 
ing Press-Schechter mass function and assuming that halos 
above a certain mass limit produce radio galaxies, and esti- 
mated that t/~6x 10" 5 . EnBlin & Kaiser (2000) however 
estimated from radio galaxy luminosity functions a total 
optical depth of the relativistic electrons to be r ~ 10~ r . 
For relativistic electrons with 7 m i„ ^ 1, one has y n th ^ T, 
and our estimate is consistent with that of Enfilin & Kaiser 
(2000). 

It is not only the cocoons of radio galaxies that have 
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been observed in X-rays, but the jets have also been re- 
cently detected by Chandra (e.g., Schwartz 2002). Although 
some uncertainties remain in identifying the emission mech- 
anism behind the X-ray radiation, inverse Compton scatter- 
ing of CMB photons is thought to be one of the possibilities 
(Stawarz et al. 2004). it is possible that the X-ray radia- 
tion from the jets may contaminate the emission from the 
cocoons and make the study of cocoons in X-rays difficult, 
especially at high redshift where angular resolution may be- 



come an issue. 



5 CONCLUSION 



We have used an evolving model of FR-II radio galaxy co- 
coons that takes into account the stopping of the jet after a 
jet lifetime of tj , and we have studied the X-ray property of 
cocoons from inverse Compton scattering of CMB photons. 
We have shown that the the X-ray power (in the 1-5 keV 
band) and the X-ray surface brightness decrease rapidly af- 
ter the stopping of the jet, as well as the radio power. At 
the same time, the X-ray photon index is found first de- 
crease to some extent while the jet is active. This occurs 
because at the important energy scale for inverse Compton 
scattering of CMB photons to keV range (7 ~ 10 3 ), two 
relevant timescales become comparable: the jet lifetime and 
that of cooling due to inverse Compton (10 7 8 yr). After the 
stopping of the jet, the photon index increases, making the 
spectrum softer. We find that typical radio galaxy param- 
eters produce an X-ray flux (1-5 keV band) of order 1-10 
/x-Crab, being emitted from regions with sizes of order of 
tens of arcsec. 

We have also shown that the ratio of X-ray to radio 
power changes little with jet power and lower cut-off in 
the electron energy distribution, and it changes mostly with 
redshift and ambient density. We have then determined the 
time-averaged ratio of X-ray to radio power and studied its 
scaling with redshift, which we found to scale as ex (1 + z) 3,8 , 
almost independent of the injected electron energy spec- 
trum, contrary to the expectations from a population of 
electrons with a single power law distribution of energy. 

Using these scaling relations, we estimated the X-ray 
luminosity function of FR-II radio galaxies at various red- 
shifts, and estimated the number of such diffuse X-ray 
sources to be ~ 25 deg~ 2 above an X-ray flux limit of 
~3x 10~ 16 erg cm~ 2 s _1 (compared to ~ 167 deg -2 dif- 
fuse soft X-ray sources observed in the Chandra Deep Field 
North for a similar flux limit). 

I am grateful to Paddy Leahy, the referee, for valuable 
comments. I also thank Mitchell Begelman for his comments 
on the paper and Biswajit Paul for discussions. 
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